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Abstract 

Normal mode phonon transmissions are studied in carbon nanotubes with the presence of Stone- 



c/3 

^ • Wales (SW) defect, using a mode-matching method and through the analysis of symmetry. The 



calculation shows that the transmission for low group velocity acoustic phonons is evidently reduced 
at high frequency range, and that this SW defect induced symmetry breaking strongly inhibits the 



'^ ■ transmission of long wave optical phonons in carbon nanotubes. The characteristic features of 

O . transmission for each phonon mode depend on the symmetry. These findings suggest that the local 



heating in the defective nanotubes may be contributed mainly from the low group velocity acoustic 



^ ■ phonons and optical phonons near the F-point. 
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Carbon nanotubes (CNTs) are uniquely symmetrioi'^'^ quasi-one-dimensional structures 
with remarkable electrical, mechanical, and thermal properties for various engineering 
applications.-i^i^i^i^ However, such symmetries for perfect carbon nanotubes are usually 
broken^, either for the purposes of engineering application such as the junction structure, 
or due to the technical flaws during the process of synthesis such as the Stone- Wales (SW) 
defects. When the symmetry of structures is broken, the equation of motion that underlines 
the physical properties of material does not respect the symmetry any more.- Therefore, such 
symmetry breaking in nanotubes is important for deep insight into both physical proper- 
ties and transport behavior. The total transmissions for all possible normal mode phonons 
at a given frequency have been first studied using the nonequilibrium Green's function 
method (NEGF).^ Of particular interest in this letter is the single- mode phonon transmis- 
sion of CNTs in the presence of the SW defect, which is known as the explicit symmetry 
breaking mechanism. Understanding single-mode phonon transmission in the symmetry 
broken nanotubes is not only fundamental to addressing the heat dissipation in nanodevices 
which challenges the current electronic miniaturization and the thermal transport in carbon 
nanotubes,-'^'^'^ but also relevant to the phonon control problem in the emerging phononic 
devices.— i^i^ This letter aims to study the behavior of single-mode phonon transmission in 
the symmetry broken nanotubes with the SW defect through the analysis of its symmetri- 
cal property. We first introduce the mode-matching method^i^ii^ used in our calculations. 
Then the phonon dispersion of the single-walled achiral carbon nanotube (11,0) is calcu- 
lated from the force constants derived from the second-generation Brenner potential.— Each 
branch of normal mode phonon is further analyzed within the full symmetry of line group 
theory. Finally we present the calculated results of single-mode phonon transmission in the 
symmetry broken CNT(11,0) in the presence of the SW defect, using the mode-matching 
method. 

The NEGF->^>^>^ has been employed to calculate phonon transmission coefficient which 
is the sum of all the normal modes from the leads. It is in principle at least applicable 
to both the elastic and nonlinear scattering regime.— However, it is difficult to analyze 
the contribution from each mode, respectively. The mode-matching methocU^'^'i^ serves 
as another equivalent way to calculate the transmission coefficient in the ballistic limit. 
In comparison with NEGF, the most prominent advantage of the mode-matching method 
is that the single mode phonon transmission can be obtained from this method explicitly. 



But the mode-matching method is vahd only in the elastic regime, where the anharmonic 
scattering is neglected. In the low-dimensional carbon nanotubes, anharmonic scattering 
does not play an important role in thermal transport below room temperature and thermal 
transport is found almost ballistic through the experiments.^'^ The dominant scattering 
mechanism is elastic due to topological defects or impurities in CNTs at low temperatures. 
Therefore, we focus our discussion on the elastic regime. The mode-matching is the exact 
method of calculating the single mode phonon transmission for the elastic scattering. 

We assume that one normal mode^'^'i^ U/„(ci;,q) = e„e**^'''^'~'^*'' is incident from the 
left lead, where each eigenvector e„ satisfies the dynamic equation of motion^ De„ = cj^e^. 
Here D is the dynamic matrix for the unit cell on the perfect lead, / denotes the atom in the 
I unit cell, R^ denotes the position for the / unit cell, and n refers to the polarized phonon 
branch. Angular frequency is denoted as uj, wavevector as q. Due to the defect scattering 
in the central region, the solution for the left/right perfect leads can be written as^ 

uf = Ui^n{uJ,q)+'^t^tui,n'{uJ,(^), (la) 

n' 

uf = 5^t^,f>,n"('^,q"), (lb) 

n" 

where n, n', n" refer to the different polarized branches of incident, reflected and transmitted 
waves. Wavevectors for the incident, reflected, and transmitted waves are q, q', and q", 
respectively. The superscript L and R indicate the left and the right. In these equations, 
t^,f„, t^4 ^^6 the amplitude transmission/reflection coefficients from mode n on the lead L to 
mode n" on lead R and to mode n' on lead L. The wave vectors q' and q" are found to satisfy 
oj = c<j„/(q') = c<j„//(q") using the generalized eigenvalue decomposition method.— Note that 
frequency does not change because the scattering is elastic. The group velocities Vg for the 
reflected waves u/^„/(ci;, q') are backward Vg < 0, for the transmitted waves Vg > 0. In order 
to obtain the coefficients t^,f^ and t^ L, s ua alternative, but equivalent, form of expression 



of Eq. ([T]) has been proposed in Ref. |13] by replacing the summation in Eq. ([T]) with the 
operation of matrix multiplication. Eq. ([1]) can be written in terms of matrix as 

uf = ii/,„(cu,q) + E(-)AK-)t', (2a) 

uf = E(+)AK+)t". (2b) 

Here E(— ) = {ei, e2, ■ ■ ■ , e„/} is the matrix formed by the column eigenvectors for the re- 
flected normal modes. A/(— ) is the diagonal matrix with the diagonal terms ranging from 



giqi'Ri^ giqjR!^ ■ ■ ■ to e*'*"'' ' . The column vector t' consists of the amphtude reflection coeffi- 
cients t^4 "with the different indices n'. Similar meaning holds for the notations E(+), A;(+) 
and t" in the right lead, respectively. With the help of Eq. ([2]), the propagation of the phonon 
waves^^'i^'i^ from / to / + s unit cell can be described by the propagator F'^ as u;+s = F'^ ■ u;, 
where F^ = E- As-E'^~^\ The inverse matrix E'^"^-' of the matrix E is pseudo-inverse because 
E is not necessarily a square matrix. 

Combining Eq. ([2]) and the propagator F^ with the equations of motion in the scattering 
region, we can solve^^ uf and uf in Eq. ([1]). Then the amplitude transmission/reflection 
coefficients t" and t' are gotten through the formulas t" = AY^{+) ■ E^^(+) ■ uf and t' = 
A'[^{—) ■ E^^(— ) ■ {ui^n — uf). The energy transmission 7^ (cj) /reflection TZni^) is related 
to the amplitude transmission/reflection^^'^ as Tn{uj) = J2n" l^n'^nP ' "^n/'/^n ^^^ TZni^) = 
J2n' l^n'nP ' '^n'/'^n) respectively, where v^ and v^ are the group velocities for the left and 
the right. 

The achiral carbon nanotube (11,0) with the SW defect is optimized with the second 
generation Brenner potential.— The force constants are derived from the optimized structure 
under small displacement. The phonon dispersion is calculated using the force constants in 
the unit cell on the lead. The calculated results are shown in the picture (A) in Fig. [H 
We further analyze the symmetrical properties for each normal mode. The achiral zig-zag 
carbon nanotube (n, 0) has the symmetry of the full line group^i^ Lz = Tl^Dnh, where T^^^ 
is the screw operation and I^nh is the point group. However, it is not necessary to work with 
the full line group. Instead, the point group is sufficient for the analysis of phonon modes 
because the normal mode vectors at the F point (g = 0) always transform as irreducible 
representations of the isogonal point group. The point group isogonal to the line group, 
i.e. with the same order of the principal rotational axis, is Dnh for achiral tubes. The 
isogonal point group for the zigzag (11,0) CNT is T)22h- Each normal mode will form a 
basis for an irreducible representation of the corresponding line group. The symmetries for 
normal modes are found by decomposing the representations of symmetrical operation into 
irreducible representations. The 12n branches of phonon modes for the zig-zag (n, 0) CNT 
can be decomposed into the following irreducible representations^ at F points 

n-l 
p«9-.a3 _ 2(Ai3 © A2u © ^Sp © ^ij © A2g © B^g © A,^ © ^s^ © J] 3{Ejg © Ej^). (3) 

We find that the calculated number of symmetries for (11, 0) CNT from the force constants 
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FIG. 1: Phonon dispersion and the symmetry of phonon modes for the achiral (11,0) CNT: (A) 
Phonon dispersion and the symmetrical notations for each mode at T points at low frequency. (Bl)- 
(B4) Vibrations of the ^2g, ^igj -^3m ^^^ E^u mode, respectively. Frequency values indicated in 
the figure are at F points. 



derived from Brenner potential strictly respect the predicted isogonal 'D22h symmetry. We 
list the point group notations at F points for the normal modes at low frequency in the 
picture (A) in Fig. [H It can be seen from Fig. [1] that symmetry properties for the four 
acoustic branches are: the longitudinal acoustic mode (LA) -^ A2U, the doubly degenerate 
transverse acoustic mode (TA) — > Ei^, and the twist acoustic mode (TW) -^ A2g. The 



calculated frequency at F point for the Aig radial breathing mode (RBM) is 254 cm ^, while 
the fitted value for RBM mode^ is ujRBM{r) = cU(io,io)(r(io,io)/r)^-°°^^^°-°°^ = 260 cm'^. Here 
^(10,10) and ?"(io,io) are, respectively, the frequency and radius of the (10, 10) armchair CNT, 
with values of ci;(io,io) = 165 cm~^ and ^(io,io) = 6.6785 A. The calculated RBM frequency 
agrees with the predicted values. The calculated vibrations for some typical normal modes 
are plotted in the picture B in Fig. [H It can be seen that the phonon modes calculated from 
Brenner potential respect the symmetries of CNTS well. 

The calculated single mode phonon transmissions for the LA, TA, and TW acoustic 
branches are shown in Fig. [2] by the solid lines. The group velocity for each mode is also 
plotted in Fig. [2] with the dashed lines. It can be seen from Fig. [2] that the transmission 
for each acoustic branch starts approximately at the value of one near the F point, which 
means that the SW defect CNTs pass long wave acoustic phonons. With the increase of 
frequency, the group velocity decreases for each acoustic mode. The transmissions for LA, 
TA, and TW modes reduce in tendency. Especially, a gap in LA mode transmission appears 
at frequency oi uj = 250 cm~^, where the corresponding group velocity decreases to a very 
small value. Compared with that of LA and TA modes, the transmission for TW mode 
has more rich features. We think that this may be ascribed to the higher symmetry of the 
atomic vibrations for the TW mode as shown in the picture (Bl) in Fig. [H 

The transmissions as a function of frequency for optical modes are shown in Fig. [31 In 
contrast to the transmission for acoustic modes, the optical mode transmissions have very 
small values at F points because of the low group velocities. This is quite different from 
the transmission for acoustic modes. The long wave optical phonons are scattered strongly 
by the SW defect in CNTs. Thus, the optical phonons at F points transport less energy 
in the SW defect CNTs. Optical mode phonons near F points contribute more to local 
heating in the SW defect CNTs. We suggest that this feature can be verified experimentally 
by measuring the intensities of Raman- active modes such as Aig. With the increase of 
frequency, the optical mode transmissions increase as shown in Fig. [3l In comparison with 
the transmission of Esu mode, the transmission for the RBM mode shows more peak features. 
This can be explained by differences in their symmetrical properties. As shown in the picture 
(B2), (B3) and (B4) in Fig. [H the vibrations of Aig mode have higher symmetry than those 
of E^u mode while the vibrations of the E^^ and E^,^ modes seem to distort the symmetrical 
structure of CNTs. Therefore, the transmission for the Aig mode is more sensitive to the SW 
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FIG. 2: Frequency dependence of phonon transmissions and group velocities for acoustic modes: 
(A) LA mode with the symmetry of A2u', (B) TA mode with the symmetry of Ei^, (C) TW mode 
with the symmetry of A2g. 



defect in CNTs, which breaks the symmetry. Together with the transmission features for the 
acoustic phonons with higher symmetry such as the TW mode, it can be inferred that the 
transmissions both for the acoustic and optical phonons with high symmetrical properties 
are sensitive to the symmetry breaking. This symmetrically mode-dependent behavior for 
the single mode phonon transmission maybe help to distinguish the different phonon modes 
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FIG. 3: Frequency dependence of phonon transmissions and group velocities for optical modes: 
(A) E^u mode with iv = 105 cm~^ at T point; (B) Aig mode with u = 254 cm~^ at T point; (C) 
E2g mode with to = 533 cm^^ at F point; (D) E^u mode with uj = 576 cm~^ at F point. 



through their transport features after the symmetry breaking. 

In summary, the single-mode phonon transmissions are calculated for the symmetry bro- 
ken CNT (11,0) in the presence of the SW defect using the mode-matching method. The 
symmetries of normal mode phonons are analyzed in comparison with the phonon transmis- 
sions. It is found that the acoustic phonon transmissions decrease with the reduced group 
velocity at high frequencies while the optical phonon transmissions have very small values 
near F points. The features of single-mode phonon transmission after symmetry breaking 
are related to the symmetry property of each normal mode. These findings may help to 
reveal the local heating problem in the defect CNTs. No evident mode conversion has been 
observed during our calculation. We also calculated other different chirality CNTs in the 
presence of the SW defect and similar findings can be made. Our present calculation holds 
at moderately low temperature where anharmonic scattering does not play an important 
role. When temperature is sufficiently high, nonlinear scattering should be considered. 
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